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Abstract 

We use mesoscale simulations to study the depinning of a receding contact line on a superhy- 
drophobic surface patterned by a regular array of posts. In order that the simulations are feasible, 
we introduce a novel geometry where a column of liquid dewets a capillary bounded by a super- 
hydrophobic plane which faces a smooth hydrophilic wall of variable contact angle. We present 
results for the dependence of the depinning angle on the shape and spacing of the posts, and dis- 
cuss the form of the meniscus at depinning. We find, in agreement with [171 ]. that the local post 
concentration is a primary factor in controlling the depinning angle, and show that the numerical 
results agree well with recent experiments. We also present two examples of metastable pinned 
configurations where the posts are partially wet. 
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I. INTRODUCTION 

Hydrophobic surfaces patterned by micron-scale posts become superhydrophobic [l|. 
Drops resting on the posts, in the Cassie or fakir state [2j, can have contact angles ap- 
proaching 180°, and roll surprisingly easily. An increasing number of actual or potential 
applications for superhydrophobic surfaces are under investigation: these include energy 
efficient condensers jsj], water harvesting devices 4], and stay dry surfaces [l|. 

Contact angle hysteresis is a term describing the difference in the advancing and receding 



contact angle of a drop as it moves across a surface [5[ . In general the hysteresis depends on 
both the drop velocity and on its shape which is, in turn, determined by how it is driven. 
Viewed as surfaces with well controlled roughness, superhydrophobic substrates provide 
an interesting model system for investigating contact angle hysteresis. Moreover a better 
understanding of hysteresis on such surfaces may be useful in their development as devices. 
For surfaces patterned with posts the advancing contact angle is ~ 180° [6] and the fluid 

PR 

motion is belived to be controlled by the receding contact line [3, [8J. Therefore in this article 
we present numerical results describing the motion of a contact line as it recedes across a 
superhydrophobic surface. In particular we are interested in how the geometry of the posts 
determines the receding contact angle at which the drop starts to move. 

Several authors have discussed how different lattice geometries and post shapes can af- 
fect hysteresis on superhydrophobic surfaces The importance of the deformation of 
the triple line , [18( has been highlighted, although a full understanding is so far 
lacking. There are recent experimental results exploring how a contact line depins from a 
superhydrophobic surface jcllsj]. 

We aim to present simulations of contact line depinning to compare to the experimental 
results and theoretical approaches. The fluid is described by a Ginzburg-Landau free energy 
model and the hydrodynamics by the Navier-Stokes equations. Hence our work is relevant 
to length scales above ~ 50nm where fluctuations are not dominant and, because we neglect 
gravity, drops smaller than the capillary length ~ 1mm. Such mesoscale simulations of 
three-dimensional drops on posts are very demanding and we found that it was not possible 
to use enough posts beneath the drop to obtain satisfactory results because of unphysical 



dependences on the initial position of the drop and on inertial effects. Therefore we describe 
a novel capillary dewetting geometry, where one wall of a microchannel is superhydropobic. 
The facing wall is smooth, and has a variable contact angle which can be tuned to drive 
depinning. 

This is described in the next section, together with details of the model fluid that we 
will consider. Then, in Sec. IIII Al we present results for the depinning angle on several post 
geometries and, in Sec. IIII B| interpret these in terms of the interface shape as it depins. 
Choi et al. [17| have recently proposed a phenomenological description of the dependence 
of depinning angle on geometry. In Sec. IIII CI we show that this is usually a very good 
representation of the data and we point out when, and why, it can fail. In Sec. IIII Dl we 
show that the simulations agree well with experiment. Finally, in Sec. |IV]we summarise and 
discuss our results. 



II. MODEL AND METHOD 
A. Simulation Geometry 

The system we investigate is shown in Fig. [TJ A channel is formed by two planes lying 
perpendicular to the z-axis at z = and z = L z . The channel is filled with coexisting 
liquid and gas such that an interface between the two phases lies within the channel. In 
the simulations the interface moves in the x-direction and periodic boundary conditions are 
imposed along y. 

Each end of the channel terminates in a reservoir of fluid and we take periodic boundary 
conditions along x (such that a second interface lies within the reservoir). This equalises the 
pressure of the gas reservoir and the gas in the channel. In the reservoir periodic boundary 
conditions are applied along both y and z. This produces a flat meniscus which does not 
exert a Laplace pressure on the fluid in the channel. A similar geometry (but with chemical 
pinning) has recently been used as a valuable tool to measure the equilibrium contact angle 
in molecular dynamics simulations J2I]. 

The bottom surface of the channel is superhydrophobic. This is achieved by patterning 
the plane by posts, of height h — 12 and rectangular cross section (d x , d y ), and choosing 



a hydrophobic equilibrium contact angle 9 eq = 100° 



, |8j. The spacing between the posts 
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FIG. 1. (a) Side and top view of the simulation geometry. A superhydrophobic surface patterned 
by square posts with 9 eq = 100° [s] faces a hydrophilic smooth channel with contact angle e q,wall- 
A slab of liquid lies between the top of the posts and the wall opposite and is terminated in a liquid 
reservoir (on the right in the Figure). A solid ridge separates the channel from the reservoir. We 
use periodic boundary conditions in both the x and y directions and, in the reservoir, also in the 
z direction, (b) Three dimensional view of the channel. 



is s x — d x and s y — d y along x and y respectively so that the fraction of surface covered 
by the posts is = d x d y /(s x s y ). We choose s x = s y = 60 and vary d x and d y to measure 
the effect of post concentration on the depinning. Between the superhydrophobic substrate 
and the reservoir we place a ridge of the same height as the posts, h, but which extends 
across the channel (see Fig. [T]). This avoids the water in the reservoir wetting the base of 
the superhydrophobic surface (to give a Wenzel state 22^ ) . 

The surface at z = L z is chosen to be smooth and hydrophilic, with a variable equilibrium 
contact angle # e q,waii- We take L z between 50 and 80 to ensure that the liquid-gas interface 
is planar in the vicinity of this wall. For high values of ^eq^aii the liquid dewets the channel. 
We define the apparent contact angle that the interface makes with the superhydrophobic 
surface as 



7 w 



180° — 6> eqiWall (1) 



(see Fig.[T]). Note that 8w is measured sufficiently far from the triple line that the liquid-gas 
profile has become flat. The receding angle Or is then the value of 0w at which depinning 
takes place. 

The size of the simulation box along z is L z . Along y the periodic boundary conditions 
allow a box of width L y = s y . At high post concentration <p it is sufficient to consider 
two rows of posts along x, such that the front is pinned to the second, in the direction 
of decreasing x. This gives a simulation box size along x, L x ~ 200. At smaller post 
concentrations we observed that Or could be affected by the bounding ridge. Therefore a 
longer system was needed with three rows of posts, as shown in Fig. [H corresponding to 
L x ~ 300. 

We turned to using the channel geometry to investigate the receding contact line after 
attempting to perform simulations with cylindrical drops lying on square posts and driven 
by a body force. These proved very difficult because limitations in system size meant that 
only a few posts could be used in the direction of motion of the drops. Inertial effects as 
the drop depinned led to results strongly dependent on the initial position of the drop with 
respect to the underlying posts. 



B. The model and numerical details 



We use a binary fluid model (with components A and B) described by a Landau free 
energy ^ written in terms of the total density distribution n = + rig and the order 



parameter ip = — n B 



23, 



* = J{ip b + \ (^) 2 ) dV- J h-<pdS (2) 

where V is the volume and S the surface of the computational domain. The bulk contribution 
to the free energy density 

^ = ynlnn + A(-^ 2 + ^ 4 ) (3) 

includes an ideal term in n which controls the compressibility of the fluid and a quartic 
contribution in ip which leads to two coexisting bulk phases with (p — ±1. The gradient 
term in Eq. (E]) accounts for the excess free energy associated with an interface. It is related 
to the surface tension between the two phases by 7 = ^8kA/9 and to the interface width 
through f = [2J. We use A = 0.02 and k = 0.04 which give f = 1.41 and 7 = 0.0267. 

The surface term controls the equilibrium contact angle 8 eq of the fluid on the solid substrate 
25| 



h = V 2kA ■ sign ( ^ - 6 eq I ■ -i /cos ( ^ ) l-cos(^) , (4) 
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where a = cos 1 (sin 2 ^ eq ). The functional (j2J) is discretised on a cubic lattice with lattice 
spacing Ax, and At is the simulation time step. The lattice velocity is then defined as 
c = Ax /At. 

The hydrodynamics of the binary fluid is given by the continuity equation (jSJ), the Navier- 
Stokes equation (jH]) and the convection-diffusion equation (J7J) 

d t n + d a (nv a ) = , (5) 
<9 t (ra; Q ) + dpinVaVp) = -d p P aP + ^[^(d^ + 9 a u^ + 5 Q/3 9 7 w 7 )] , (6) 
d t <p + d a (<pv a ) = MV 2 /i , (7) 

where we have introduced the flow velocity v and the dynamic viscosity 77. In Eq. the 
parameter M is a mobility which controls the diffusion of the triple line when the system is 
pushed out of equilibrium. This is the mechanism by which the contact line moves although 



no slip boundary conditions (v = 0) are implemented at the solid wall. The pressure tensor 
P a p in Eq. ([6]) and the chemical potential /i in Eq. (J7J) are derived from the free energy ([2]) 
and are given by 24] 



P. 



fi = A (—ip + y? 3 ) - kVV, 
.. < = (Pt - f (9 7 ¥>) 2 - Mpd^y) 5 a/3 + K(d a (p)(dp<p) 



(8) 
(9) 



The mesoscopic equations (jM7|) are solved using a lattice Boltzmann algorithm 26) in the 
implementation described in detail by Pooley et. al. [27]] which has been shown to be suc- 
cessful in reducing spurious velocities related to lattice artifacts. For further details we refer 



to Refs. 
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28) . Although our results have been obtained using a binary model, we 



will be working in a quasi-static regime where a liquid-gas formalism gives the same results. 
We have chosen to use this approach because our experience is that it relaxes efficiently to 
stable and metastable minima of the free energy and because it can very naturally be ex- 
pended to dynamical problems when the drop is moving across the surface. Surface Evolver 



29) provides an alternative approach in the static limit 30M32 ] . 



III. RESULTS 



Fig. [2] summarises the simulation results for the depinning angle Or for several different 
post geometries. After discussing the data, we relate it to the shape of the meniscus as the 
interface is pulled across the posts. We then interpret the results in terms of a phenomeno- 
logical model of depinning, and compare to experiments on drops on superhydrophobic 
surfaces. 



A. Post geometry 

Ridges across the channel. Before considering a substrate patterned with a more 
complex post geometry it is instructive to sketch out the bounding cases of ridges and 
grooves either running across (s y = d y in Fig. [TJ or along the channel (s x = d x ). In the 
first case the system is translationally invariant along y and, noting that the interface is 



flat, Or follows easily from the Gibbs' criterion 



33) . This states that dewetting occurs if 




FIG. 2. Dependence of the receding contact angle Or on concentration. We consider three 
geometries: square posts arranged on a square lattice (x), rectangular posts with d y = 2d x on a 
square lattice (□), and square posts on a square centred lattice (A). In each case Or lies between 
the two values shown. For the larger /smaller of these the front dewets/remains pinned. The light 



lines are the scaling predictions of Ref. 



171 ] . The heavy lines refer to ridges and grooves along x 



(Or = 0c, full line) and y (Or = 6> eq , dashed-dotted line). 

the angle between the meniscus and the wet, horizontal face of the ridge is smaller than 
the equilibrium contact angle of the flat surface. Hence Or = eq independent of <fi. The 
simulations reproduce this limit, shown by the dash-dotted line in Fig. [2j 

Ridges along the channel. The case of grooves and ridges along x, parallel to the flow 
direction, is slightly more involved because the liquid-gas interface is not planar. In this case 
there is no pinning and the channel is spontaneously wet or dewet for Ow > Or or Oy/ < Or 
respectively, while for Ow = Or the liquid-gas interface remains stationary. The shape of 
the interface is independent of its position in the channel and therefore, withdrawing the 
interface through 5x, the change in free energy is 

5 J 7 oc —7 • Sx ■ cosOw + 7 • Sx ■ ((f> ■ COS^cq + — 1). (10) 



At the threshold angle 5 J 7 = 0. Hence Or = 0c where cos 0c = 0cos6* cq + 



[2 



and 0c can 



be identified as the Cassie angle of the superhydrophobic plate as expected [2|. This limit 
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is shown by the full line in Fig. [2j 

Square posts. We now consider a bottom plate patterned with square posts. We vary 
the post size, d x = d y while keeping the period of the pattern, s x and s y , fixed. Fig. [2] shows 
numerical estimates of the receding angle Or for different post concentrations (p. For a given 
(f) we ran simulations at different Ow- In Fig. [2] we bound Or with the two values of 6w 
between which the crossover from pinning to dewetting was observed. The results show that 
Or tends to eq for — > 1, and to 180° with a divergent slope for — > 0. Thus the numerical 
results for Or are, as expected, bounded by the depinning lines of the grooved geometries. 

Rectangular posts. In Fig. [2J symbols (□) indicate simulation results for rectangular 
posts with d y = 2 ■ d x arranged on a square lattice. As expected depinning is harder than 
for square posts for a given concentration, primarily because the spacing between posts 
along y is smaller and the geometry interpolates between squares and ridges along y. As 
<p — > 0.5, Or — > eq , and continuous ridges form across the channel. Similarly, rectangular 
posts elongated along the x direction will make depinning easier at a given concentration. 

Face-centred square patterning. Finally we present results for a different arrange- 
ment of posts, square posts arranged on a face-centred square lattice. This geometry is 
shown in Fig. [2j where the triangles depict the corresponding simulation data for the de- 
pinning angle. For a given <p it is easier for posts on a face-centred square lattice to depin 
than posts in a square array because of the larger distance between the posts. However, 
for a given post spacing, the opposite holds true because the post in the middle of the unit 
cell helps to pin the receding interface. This effect, which becomes dominant for the largest 
concentration reported in Fig. [21 will be discussed in more detail in Sec IHI Bl 

B. Menisci deformations 

In general on patterned surfaces the receding contact angle is larger than the Young 
angle because of the deformation of the liquid-gas interface near the triple line. As the 
liquid bulges between two posts, as shown in Fig. [3], the interface in the vicinity of the post 
corners is pulled forwards thus decreasing the local contact angle and facilitating depinning. 
This deformation is strongly affected by details of the patterning geometry. It occurs in the 
vicinity of the posts, and the channel height L z is chosen to be sufficiently large that the 
interface is essentially planar near the top of the channel. This ensures that the depinning 
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FIG. 3. Pinned menisci recorded one liquid layer above the posts in the xy plane at (a) constant 
post concentration cj) for different values of the apparent contact angle Ow and (b) at Ow ~ Or 
for different values of 0. (c) Pinned menisci in the xz plane, comparing cross sections which run 
through the centre of the posts, and half-way between the posts for different Ow and <ft = 0.17. (The 
pinned interfaces do not appear to lie exactly at the corners of the posts because we record their 
position one lattice spacing above the surface and because of the diffuse nature of the interface.) 
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angles are independent of L z [34]. 



In Fig. [3(a) and (b) we show the profiles of pinned menisci in the xy plane one lattice 
spacing above the top of the posts. The first of these plots shows that, as Or is decreased, the 
meniscus remains pinned on the back of the posts, but is gradually pulled forwards between 
the posts. Depinning occurs at Or = 125°. (In the Figure the pinned interface does not 
appear to lie exactly at the post edges because the interface is diffuse.) 

Fig. [3(b) shows how the deformation between the posts depends on different values of 
the concentration </> for $w ~ Or. As <fi decreases the larger gaps mean that it is easier 



for the interface to deform, leading to a 



ligher Or. Moreover the interface shape become 



7, 



increasingly cos /i- like, in agreement with 

In Fig. 0(c) we compare the shape of the menicus in the xz plane, taking two cross 
sections, one cutting the centre of a post, and one lying halfway between two rows of posts. 
As 8\v — > Or the difference between the two profiles increases, meaning a bigger distortion 
in the vicinity of the posts in agreement with Fig. [3(b). Notice that, because the interface in 
the reservoir is flat (Fig. [T]), the Laplace pressure is zero and the pinned liquid-gas meniscus 
has zero mean curvature. (This is consistent with drop hysteresis experiments [8( where the 
radius of the drop is much bigger than the length scale of the surface patterning.) Therefore, 
as the liquid bulges forward, the xy profile of the interface between the posts becomes convex 
(to the liquid) and hence, to preserve the condition of overall zero curvature, the xy profile 
in the centre of the posts must become concave, as shown in Fig. [3(a). As this curvature 
increases the contact angle with the top of the post decreases until it reaches eq when 
depinning can occur. 

As the post concentration is decreased the meniscus underneath the liquid does not remain 
flat but increasingly bulges between the first and second row of posts as shown in Fig. [H 
This occurs because the meniscus is curving sharply upwards to meet the top surface of 
the channel unless it is pinned by the posts. Hence upwards curvature between the posts 
contributes to lowering the interfacial free energy. Recent experiments have shown that small 



drops can be left behind on the posts as the interface moves forward 35|] and, although we 
were unable to access a parameter range where this occurred in the simulations, the distortion 
of the interface under the drop may be a precursor of such a pinch off. 

We now consider the face-centred square arrangement of posts. At low post concentra- 
tions, ((p = 0.109 and <fi = 0.22 in Fig. [2]), depinning proceeds as for square posts: once the 
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FIG. 4. Meniscus profile in the xz plane showing that the curvature under the liquid between the 
first and second rows of posts increases with decreasing concentration: (a) <f) = 0.25 (b) <p = 0.0625. 
The broken line shows the interface profile halfway between the posts and the full line the profile 
in the centre of the posts. (The pinned interfaces do not appear to lie exactly at the top corners 
of the posts because of the diffuse nature of the interface). 



interface has completely depinned from the rear edge of a row of posts the channel dewets. 
For <p = 0.347, however, we find that depinning occurs through a two steps process. This is 
illustrated in Fig. [5] where we plot pinned profiles recorded in the xy plane and in the xz 
plane. The profiles in the Figure indicate that for 6w = 132° the front remains pinned to 
the rear edge of a row of posts, as for simple square patterning. The only difference is that 
the front between two posts is also pinned to the post in the middle of the unit cell (labelled 
B in Fig. However for 128° < 6\y < 132° the contact line, instead of dewetting, remains 
pinned to the rear of the second row of posts (B), leaving the first row of posts (A) only 
partially wet. Finally for 8w < 128° the liquid completely dewets the A posts. The channel 
will then empty as 6w is lower than the value needed to drive either transition. Note that 
in Fig. [2] we give the value of Or corresponding to the second, final depinning transition. 

The results for square posts shown in Fig. [3] suggest that the corner of the posts can be dry 
when there is still pinning on the back edges. However this result could be affected by the 
diffuse nature of the interface. To demonstrate unequivocally the existence of pinned states 
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FIG. 5. Pinned menisci shapes for the face-centred square patterning for (j> = 0.347 and three 
different Ow- (a) Profiles in the xy plane one lattice spacing above the top of the posts (along the 
dotted line in Fig. [Sb). (b) Profiles in the xz plane for two values of y (corresponding to the dotted 
lines in Fig. [5h.). For 9w = 132° the interface is pinned to the A and B posts whereas for 8w = 129° 
and 9w = 128° it is pinned to the B post, leaving the A posts partially wet. In (a) the profiles do 
not, at first sight, appear to be pinned to the post B because they are recorded one lattice node 
above the posts and the meniscus meets the dry face of B at an angle which approaches 90° as 9\y 
increases. This is apparent in frame (b). 
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FIG. 6. Time evolution of the liquid-gas interface on T-shaped posts for 8\y = 135°. The profile is 
recorded one lattice unit above the top of the posts, every At = 6 • 10 4 to 39 • 10 4 lattice Boltzmann 
steps. The starting configurations are the flat interfaces labelled by (a), (b) and (c). The interfaces 
quickly relax to curved configurations and then much more slowly translate along x to a pinned 
state in which the posts are partially wet. 



with the post tops partially wet we consider T-shaped posts. We chose as initial conditions 
three flat interfaces lying at different values of x, labelled (a) (b) and (c) in Fig. |6l and 
6\y = 135°. The interface very quickly relaxed to a curved configuation which reflected the 
presence of the posts. It then advanced or receded much more slowly across the tops of the 
posts towards a final equilibrium state, lying between the full and dashed profiles in the 
Figure, that was independent of the starting configuration. The partial depinning occurs 
because 135° is bounded by the value which depins a d y /s y = 10/60 ridge {Or = 150°) and 
that which depins d y x d x = 30 x 10 rectangular posts {Or = 131°). For more complicated 
geometries, it is possible that several metastable pinned configurations could exist. 
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C. An estimate of the depinning threshold 



In a recent paper Choi et al. 17j proposed a phenomenological theory for the depinning 
angle of a receding contact line. They argued that the receding contact angle can be esti- 
mated using a Cassie-like argument, analogous to that described in Sec. IIII Al for a surface 
patterned by ridges and grooves along x, but using the local post concentration 0£ sampled 
by the interface as it advances a small distance beyond pinning 

cos 6 R = 4> e cos 6> cq + 0£ - 1 . (11) 

Eq. (jlip assumes that, when the interface depins, it preserves its shape as it moves forwards. 
Predictions of the scaling relation (fTTj) . for each of the simulation geometries we consider, 
are reported in Fig. [2] as light lines. 

For square posts the linear density sampled by the interface is (pi = \/0. As is apparent 
from Fig. [21 the estimate ffTTl) is in excellent agreement with the numerical results over 
a wide range of 0. Agreement is essentially exact at high post concentrations, where the 
meniscus deformation is confined to within the first line of wetted posts. At low 0, deviations 
from Eq. (TTT|) start to be apparent, in particular at the lowest post concentration simulated 
= 0.0333. This is because the meniscus is more strongly deformed and probes the region 
of the channel with no posts. Hence it changes its shape as it moves forwards. 

It was not possible to obtain simulation data for lower because this would require 
too big a simulation box. However, the discrepancies might be expected to be larger for 
rectangular posts elongated in the y direction which, for a given 0, extend less far along x. 
Choosing, as in Fig. [21 d y = 2 • d x gives 4>g = \^2cf>. The agreement is still pleasing, but less 
close than for square posts. Eq. (jTTj) gives a pinning angle which is slightly too high, and the 
increase in the discrepancy with decreasing is clearly seen. This occurs because depinning 
is facilitated by the detailed shape of the interface between subsequent rows of posts, which 
is not accounted for the theory. 

By contrast, for the face-centred square post geometry the prediction of Eq. ( II ip (taking 



<pe = a/ 0/2) underestimates the depinning angle, and becomes much less accurate for large 
0. This is because, as argued in Sec. IIII BL the post in the middle of the unit cell helps to 
pin the receding interface. In particular, for the highest post concentration = 0.374, Eq. 
(fTTj) predicts a 9r which is compatible with the first dewetting transition of Fig. [51 rather 
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FIG. 7. Contact angle hysteresis as a function of the post concentration eft. Circles: experimental 
results from [g]. Triangles (digitised): experimental results from [8|. Crosses: simulations for a 
surface patterned with square posts. Broken line: scaling prediction from [17]. Experimental errors 
are of order the scatter in the data. 



than reproducing the second depinning which finally dewets the channel. 



D. Comparison to experiment 



We next compare the simulation results to experiments measuring drop hysteresis on 
superhydrophobic sur: 



Reyssat and Quere 



,8 



nevs 0,0]. In two independent sets of experiments, carried out by 
and by Priest at al. drops were placed on a superhydrophobic 
surface which was then tilted. The advancing and receding contact angles, 9 a and 9r, were 
measured at the point where the drop started to move, and the data for each experiment 
was reported by plotting A cos 9 = cos 9a — cos 9r as a function of the post concentration. 
In the experiments described in Ref. 8| the superhydrophobic surface was patterned with 
circular posts of diameter 2.4 fim and had an equilibrium contact angle 100°. The radii of 
the drops was between 1.3 and 1.7 mm. In Ref. (|, by contrast, the posts were square with 
sides of 20 /im, the surfaces had an equilibrium contact angle between 95° and 115° and the 
drop radii were < 2 mm. 
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There are differences between the experimental protocol and our simulations. Firstly 
we consider a slab of liquid which is translationally invariant along y sufficiently far from 
the posts, rather than a drop which forms a spherical cap. Secondly in the experiments 
the drop is pushed by gravity, a body force, rather than by a capillary force due to the 
changing contact angle of the upper wall of the microchannel. These should not be important 
differences for drops of the size considered in the experiments. Thirdly we assume in the 
simulations that depinning is controlled by the receding contact angle, and we do not measure 
an advancing contact angle. For the purposes of comparison we set 6 a = 180°. Other sensible 
choices (eg. in jg] 9a is measured to lie between 160° and 170°) make little difference to the 
estimate of A cos 9. 

Fig. [7] shows that simulation results compare well with experiments, in particular for the 
square posts |6|. Indeed the agreement is surprisingly good given the differences between the 
experimental and simulation geometries. This supports the assertion that drop hysteresis 
on superhydrophobic surfaces is controlled by the receding contact line 



IV. CONCLUSION 

We have introduced a novel simulation geometry, a microchannel with a superhydropho- 
bic wall. Considering capillary dewetting of the channel has allowed us to obtain numerical 
results for the receding contact angle 9r of several superhydrophobic surfaces. We used suf- 
ficiently high channels that 9r only depends on the details of the superhydrophobic surface, 
and therefore we expect our results to be relevant to drop geometries, as long as the drops 
are sufficiently large, and as long as depinning is controlled by the receding contact line. 
Comparison to experiment indicates that this is indeed the case. 

Our main results are summarised in Fig. [2] which shows the receding contact angle at 
depinning as a function of concentration for square and rectangular posts arranged on a 



simple square or a face-centred square lattice. We find, as first discussed in [17] . that the 
primary control parameter for depinning is d y /s y , the fraction of space occupied by the 
posts in the direction perpendicular to the forcing. This is because depinning is easier if the 
interface can more easily deform between the posts. Deviations from this scaling manifest 
themselves once the interface penetrates beyond the posts, for example, for square posts 
at low concentrations, or if it reaches a second row of posts, as for face-centred square 
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patterning. 

The existence of pinned menisci with partially wet posts has been an open question 



We find that this is possible for posts arranged in a face-centred square lattice where the 
interface depins in two steps. We have also identified a metastable pinned state for which 
T-shaped posts are partially wet. However, detecting pinned menisci where the posts are 
partially wet in experiments is likely to be challenging due to inertial effects. 

From a methodological standpoint, we hope that the geometry we consider here will prove 
useful in future simulations of depinning. Open questions which we will aim to address 
include dynamic contact angle hysteresis - measured when the liquid is moving - and the 
response of a moving slab of fluid to a driving force provided by changing the contact angle 
of the wall. It would be of interest to fabricate an experimental realisation of the patterned 
channel used in this work, with the contact angle of the hydrophilic flat wall # e q )W aii controlled 
by an electrowetting potential. Looking to applications, contact line pinning by posts is being 
used to control fluid m o«o n in uflcrochannels fl and au unde.taudiu, of pinning may help 
design better slippery surfaces, or microchannels with anisotropic filling properties. 
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